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^ ■ Abstract 

! The polarised singlet coupled DGLAP equations (LO) are trans- 

■ formed by a Taylor series expansion at low x into a pair of partial 

differential equations in x and t(t = lnQ 2 /K 2 ). The pair of coupled 
Qh| partial differential equations is then reduced to canonical form and 

p ' the resultant system solved by applying the method of characteristics 

under small x and t approximations. The result is compared with 

. . ! some exact solutions available in the literature. 

> ■ 

b : 1 Introduction 

The measurement of the polarized structure functions of the proton by the 
European Muon Collaboration in 1988 revealed that the structure of the pro- 
ton is more profound to be understood only from the unpolarized DIS. Their 
result P] E], which is often referred as the 'proton spin problem' indicated 
that the contribution of the quarks to the nucleon spin is very small. The 
rest of the spin then must be carried by the gluon and/or by the angular 
momentum of quarks and gluons. The spin distributions among the different 
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contributors is expressed by a sum rule [3J Ej 



\= l -XL + ^G + L q + L 



(1) 



where AE is the contribution from the quarks, AG from the gluons and L q 
, L g are the orbital angular momentum of quarks and gluons to the total 
spin of the nucleon. In the static quark model with SU(6) symmetry, the 
spin is entirely carried by the quarks i.e AS = 1. In the naive quark-parton 
model, AS is predicted to be ~ 0.7. These are in contradiction to the EMC 
result. Based on the next-to-leading order QCD analysis of the deep inelastic 
polarized data, there are several model parametrization El HI IH1 HH CD] of 
the polarized parton distribution functions. All of these models predict that 
the value of AS is around 0.3 or less; that is, the quark contribution to the 
nucleon spin is very small contrary to our expectation. In view of this, the 
polarized gluon distribution has attracted in recent years special interest in 
polarized DIS, because, a rather large contribution of the first moment of 
the gluon distribution is essential to compensate for the low value of the 
quark contribution to the total helicity of the nucleon. Recent data JT] 
and analyses 13 El El ED] suggest that the contribution of the polarized 
gluon momentum is quite large, though the error involved is too large to 
significantly constrain its contribution. It is therefore imperative to know 
how the gluon is polarized inside the nucleon. In this paper we investigate 
the polarized gluon and the quark singlet contribution to the nucleon spin 
with the help of the DGLAP evolution equations. We adopt the method of 
characteristics to solve the coupled polarized DGLAP evolution equations. 
In 321 we discuss the formalism and in §21 we discuss our results. 

2 Formalism 

For Q 2 evolution it is convenient to introduce a flavour singlet combination 
of polarized quark and anti-quarks by 



where the sum is over all quark and anti-quark flavours. Here Aq(x, Q 2 ) = 
(q + (x,Q 2 ) — q~(x,Q 2 )), where q + (x, Q 2 ) ,( q~(x, Q 2 )) refers to the quark 



AS(x, Q 2 ) = J2(Aq(x, Q 2 ) + Aq(x, Q 2 )), 



(2) 
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distribution with helicity parallel (anti-parallel) to the parent proton. The 
polarized gluon distribution is similarly defined as 



Ag(x,Q' 2 )=g^(x,Q 2 )-g^(x,Q 



H, 



(3) 



The Q 2 evolution of the polarized singlet and the gluon distribution mix 
non-trivially and is given by the DGLAP ^21 Ell E] equations 



d_ 

di 



AE(x,t) 
Ag (x,t) 



a s (t) ( AP C 



2n f AP gg 



2tt \AP gq AP gg 



AE (x,t) 
Ag (x,t) 



(4) 



where t = Irtj^ and APij are polarized splitting functions which are known 
at LO HU and NLO [13113 01! • The symbol ® stands for the usual Mellin 
convolution in the first variable defined as 



(AP x Aq)(x,t) = f 1 ^AP (-) Aq(y,t). 



(5) 



The polarized parton distributions at an initial scale enter as boundary val- 
ues in the solution of the above equations. However, with exact splitting 
functions, analytic solutions in the entire range of x is not possible. The 
splitting functions can be simplified by taking their moments and expanding 
around the rightmost singularity at iV = where N is the moment variable. 
In this procedure the splitting functions[9j at LO get simplified to 



= It 1 + ¥( l - *)], ApW) (a; ) = 1 2n f [-l + 25(1 - x)}, 
a d(o)/^\ _ 4ro xn ™M a r>(0)r„\ — qIa 13 xn „m n f xn 



AP(°)(x 



(6) 



OT W 3 L- — 6 6{1-X)]-^6{1-X) j 

With these simplified splitting functions, Eqs.Q) can be solved analytically 
with some approximations about small x behaviour of these to be discussed 
later. The two equations Eqs.(jU can be written separately in LO as 



0_ 

Of 



AS(x,t) 



a s (t) r 



2tt 



AP$ ® AS(x, t) + 2n f AP$(x) ® Ag{x, t)] (7) 



and 



|A 9 (x, t ) 



a s (t) 
2tt 



AP? <g> AS(x, t) + APjJ(x) <g> A 5 (x, t) 



(8) 



Using the convolution integral(EqlHJ), we write Eq.0 and Eq.(jHJ) as 



2tt 



djiAPj^ A£(~,t) + 2n, (feAPjW Ag{-,t) 

J x Z J x Z 

(9) 
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-Ag(x,t) = -±+ 
at 2tt 



f 1 dz APW (z) A£(-,t) + f dz AP^(z) Ag(-, t) 

J x Z J x Z 



(10) 

x 

Introducing a variable u as u = 1 — z, we can write — as a series 

z 

XX ^ k 



5>- (ii) 



— = = x > ii 

2 1 — it r - ?, 

«=o 

Now we expand A£(— , t) and Ag( — , t) in Eq.tjUJ) and Eq. ffTUj) in Taylor series 
using Eq. ffTT]) as 

AE(-, t) = AE(x + x y u\ t) = A£(x, t) + (if ^ ^M) + ... ( 12 ) 
and 

(-, f) = A^Or + x £ f) = A^x, f) + (x £ u *) 9 ^ x ^ + (13) 

z fe=i fe=i ^ x 

At low a; the above series Eq. (jl!2|) and Eq. (jl3j) are convergent |18j. Hence, 
neglecting higher order terms, the series can be approximated as 

AE(* , t) « A£(x, t) + (x f; u ») g^!fc*) (14) 
z fc=l 

and 

A^(-, t) « A^x, t) + (x £ (15) 
2 fe=i ax 

Using Eq. ffnj) and Eq.fJTSJ) and the splitting functions given by Eqs.tjHJ), we 
can write Eq.© and Eq. (fT0|) as 

f^^=/('Hx)A^x, t ) + /f^|M +41 ) WA9(l , t ) +/f ^i) 

(16) 

and 

(17) 
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where we have used a s (t) 



47T 



in LO. The quantities = q,g;j 



1, 2, 3, 4) in Eq.(J]U) and Eq.flE} are given by 



/ft 



= 2n,- / _ (-1 + 2,5(1-*)) , 
4 2) (^) = 2n/i [ ^ (-1 + 2 5(1 - *)) ^ f> fc ) 



4^ 



4 3) ( 



4 r 1 dz , 



4 4) ( 



and 



/J 4) ( 



1 

x Z 



1 dz 



13 



3(4-— 6(1 -z) ) -^5(1-2 
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3(4-^(l- 2 ))-^(l - : 



\ k=l 



(18) 

(19) 

(20) 
(21) 

(22) 
(23) 

(24) 
(25) 



Carrying out the integrations in Eqs. (jl8H25j) . we recast Eq. ([TK|) and Eq. (jT7jl 
as two coupled partial differential equations in two variables x and t as : 



, dAE(x,t) , 0Ag(x,t) dAZ(x,t) , dAg{x,t) 

°11 £U + a l2 a, + °U Q „ + °12 



Of 



Of 



Ox Ox 
R' n AZ(x,t)+R' l2 Ag(x,t) 



(26) 



and 



, d£E(x,t) , 0Ag(x, t) , dAE(x, t) , 0Ag(x,t) 

a 21 oT I" a 22 oT I" "21 q I" "22" 



Of 



Ot M Ox 11 Ox 

= R' 21 AE(x,t)+R' 22 Ag(x,t), 



(27) 
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where 



a' u — t, a[ 2 — 1 /nn\ 

a' 21 = 0, a> 22 = t f ' ^ 



ft 



21 



-§(!-* 



xln- 

X 

(l -2x-x/n±) 



y 12 



-2n 



/ 2 



•1 + x + xZn- 



b' 22 = -12 (l -x-xZn±) 



(29) 



and 



n 



3 V2 



*J 1 T- 



12 



2n/| 



^2i = |(2^-l 



22 



-ln\ 
12 ini-(f + ¥ 



(30) 



Eq.(j2Bj) and Eq. (}2Tj) represent a system of first order coupled partial differ- 
ential equations. We can reduce these equations to canonical form and then 
solve by the method of characteristics. To do that, we introduce a vector 
Au(x, t) defined by 

' AE(x,t) \ 



Au(x, t) 



Ag(x,t) 



and express the two equations Eq. (j26|) and Eq.(|27|l in matrix form 

a'Au t (x, t) + b'Au x (x, t) = R'Au(x, t), 
where the matrices a' , b' and R' are 



l 21 



*12 
'22 



and 



b' 



R' 



h' b' 

u 21 u 22 



R'u R l2 

-^21 -^22 



(31) 
(32) 

(33) 
(34) 

(35) 



the elements being given by Eqs. 

non-singular, we multiply Eq. (j32j) from left by a' -1 and get 
Au t (x, t) + A' Au x (x, t) = B' Au(x, t) , 



respectively. The matrix a 1 being 

(36) 
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where the new matrices A' and B' are: 

A' = a'^b' (37) 

and 

B' = a'^R'. (38) 

Eq. (j32j) represents a system of two coupled first order partial differential 
equation in the two variables x and t for the vector Am prescribed by Eq.(|3*T|). 
Its principal part, i.e. Au t (x, t) + A'Au x (x, t) is completely characterized by 
the coefficient matrix A'. Since the matrix A' has n [here n=2] distinct 
eigenvalues, the system Eq. fJHol ) is a hyperbolic one and it is possible to 
obtain its canonical form in the following way [TTfl I2U| : Let X'^(i = 1,2) be 
the eigenvalues of the matrix A' (EqJ37|) and P' be the matrix formed by 
the corresponding eigenvectors. Now, if A is the diagonal matrix with the 
eigenvalues Ai and A2 as the two elements, then we have 

Let Av(x,t) be a new function defined in terms of Au(x,t) by the equation 

Au(x,t)=P'.Av(x,t). (40) 

so that 

Av = P'~\Au. (41) 
Differentiating Eq.pO)) with respect to t and x respectively we get 

Au t = P'Av t + P' t Av, Au x = P'Av x + P' x Av. (42) 

Substituting Eqs.flHJ) and (02)) in Eq.flUS), we obtain 

P'Av t + P' t Av + A'P'Av x + A'P' x Av = B'P'Av. (43) 

Multiplying Eq. (j43|) from left by P'^ 1 and using Eq. (j3*9*j) . we obtain 

Av t + AAv x = Ae, (44) 

where 

Ae = P'-\B'.P' - P[ - A'.P' x ).Av . (45) 
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Eg . (J44j) is in canonical form. In components, it is 



AvP + A' (i) Av« = Ae (i) (i = 1, 2), (46) 

dAv® ■ dAv® 

Eg .(14611 shows that the principal part viz. — h A'^ — involves only 

at ox 

the component Av^' of the vector Av and its derivatives i.e. the eguations are 
decoupled. These eguations (Eg l4T)J) can be reduced to ordinary differential 
eguations 

r /A ?; («) 

— ^— = Ae {i \x,t,Av {i) {x,t)) (i = 1, 2) (47) 
along the characteristic curves defined by the eguations 

V»(«,*) = ^. (48) 

In order to integrate the ordinary differential eguations (EqJ47jl to get an- 
alytical forms for Av, one has to express the right hand side in terms of t. 
This is done from the solution of the characteristic eguations Eg .([48)1. which 
expresses as a function of t. However, the integration of Eq.([48|) to get 
such an expression depends on the nature of the eigenvalues X'^ of the ma- 
trix A'. For simple forms of the eigenvalues, it is possible to integrate Eg. (J47|) 
analytically. We discuss this below. 
The eigenvalues of the matrix A' are 

X>^ = ±(b> 11+ b> 22 ±s), (49) 
where 



b' 2 1 + 4b' 12 b 21 -2b' u V 22 + b> 2 22 . (50) 



Taking the elements of the matrices a' and b' as given by Egs. (j28H29j) . the 
eigenvalues are found to be 



X /( - 1,2 \x, t) 



27t 



-10 + lOx + 10xln- 

x 



1 l 2 



±\ 55 - 101 x + 46 x 2 + (-1 10 x + 101 x 2 ) In- + 55 x 2 (In 

V xx 



(51) 
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where the + (— ) sign corresponds to A'W(A'^). And the two components of 
the vector Ae are 



+ t 22 : F'')fi2i 



9.< 



+ =F- 



(-6'u+6 22 -«) 



" 

2 s 



96' 



2fc' 21 fl' 12 (-b' 11+ b 22 Ts)R 22 
' 9s 9s 



Av 2 



(52) 



where the upper sign (-) corresponds to Ae^ and the lower sign (+) to 
Ae®. We note that both A /(i) and Ae (i) are factorisable in x and t, the 
quantities b\ • and being functions of x only. As we have already observed, 
integration of Eq. ()47|) depends on the analytical solution of the characteristic 
equation . But with the above eigenvalues (Eq.(|5Tjl). analytical solution of 
the characteristic equations(Eq. ()48Jl ) cannot be found and so we cannot find 
the necessary transformation equation to express x as a function of t for the 
integration of Eq.(|4"7)). But under certain extreme situations, it is possible 
to get analytical solutions which we discuss below. 



2.1 Solution when x —> 

The eigen values of the matrix A' can be simplified by making certain ap- 
proximations about the elements of the matrix b'. Before we do that, we 
write the eigenvalues in simple form. The eigenvalues being factorisable in x 
and t can be written as 



A'«(z,t) = \' {ia \x) X' (ib \t) (i = 1,2), 



(53) 



where 

\ ,{ia Hx) 



2 

27 



-10+ lOx + lOxln- 



x 



±\/55 - 101 x + 46 x 2 + (-110 x + 101 x 2 ) In- + 55 x 2 In- 
V xx 



(54) 



and 



v 1 t 



(55) 
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In Eg . (|54j) . +(-) sign corresponds to i=l (2) respectively. Now in the limit 
x — > 0, Eq.([54|) takes the values 

A'< fa > _> 1 (_io ± V55) . (56) 
The characteristic equation Eq. (J48j) now takes the form 

=*"•»<«)*■ w 

Integrating Eq. fjBTj) we get 

x {l \t) = a ii) lnt + C\ (58) 

where 

a {€) = — (-10 ± i = l,2 (59) 

and C« are two constants of integration. To get the constants of integration, 
let (x, t) be a fixed point ^J] i n the (x, t) plane through which the two 
characteristic curves Eq. (|58p pass through. That is, x (!)(*) =iand x (2) (t) = 
x. Then from the Eqs.()58p we get 

x {i) {t) =x + a {i) ln(^j . (i = 1, 2) (60) 

These are the two characteristic curves that pass through a common point 
(x, t) up to which we can evolve the two functions Av^ 1 ' and Av ^ defined 
by Eq. (|41|) . Furthermore, if the two characteristic curves cut the initial line 

t — to (to = at ^^(*o) = T i an d x^(t ) = t' 2 respectively, then 

r[ = x+a w ln(*-f) (61) 

and 

T i=x + a (2 hn^ . (62) 

Hence the two characteristic equations corresponding to the two eigenvalues 
are respectively 

x m(t)= T [ + a^ln(-\ (63) 
\t J 
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and 



x { - 2 \t)=r' 2 + a { - 2) In 



to 



(64) 



In this approximation, the two components of the vector Ae defined in 
Eq.(gnj) are 



Ae« 



2970 



2255 + 284 v / 55)Aw (1) + 3(275 - 52^55) Av 



,( 2 ) 



(65) 



+ (4400 - 404V55)Aw (1) + 4(-880 + 119V55) Av {2) In— ^j- 

x^>{t) 



and 



2970 



(2255 + 284v / 55)At» (2) + 3(275 + 52 v / 55)At» 



(66) 



+ (4400 + 404 v / 55)At» (2) - 4(880 + 119\/55)Az; (1) /n^-- 



Using Eq.(|55j) and Eq.([Bl)j). we rewrite the two differential equations Eqs.(|4Tjl 
separately as 



dAv^(x,t) TA , , „ , NJ 1 ,1, 
- + B 1 (x,t)ln }-dt 



and 

where 

Ai(x,t) -- 

B x {x,t) = 
A 2 {x,t) 

and 

B 2 (x,t) 



AvW(x,t) 

dAyW(x,t) 
AvM(x,t) 



xM(tyt 



[A 2 (x,t) +B 2 {x,t)ln M/ .A ~dt. 



J xW(t) s t 



2970 
1 

2970 
1 



2255 + 284V55) + 3(275 - 52V55) 



2970 



2970 



(4400 - 404V55) + 4(-880 + 119V55) 
-(2255 + 284\/55) + 3(275 + 52^55) 



Av^y 



AvW 



(4400 + 404V55) - 4(880 + 119V55) 



Av( 2 ) 
Av ( 2 ) 



(67) 
(68) 

(69) 

(70) 
(71) 

(72) 
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Eq.(|67|) and Eq. (J68|) can now be integrated. Integrating along the character- 
istics from t = to to t = i we get 



Av {1 \x,t) = Av {1 \r[) exp 
and 

Av (2 \x,t) = At/ 2 ) (7-2) exp 



to 



A x {x{t),t) +B 1 (x(t),t)ln 



dt 



xW(t) t 



(73) 



A 2 {x{t),t) + B 2 (x(t),t)ln 



1 1 dt 

xVHfjfT 



Or changing the variables from (x, t) to (x, t) 



(74) 



Av^{x,t) = Av^{t[) 



exp 



to 



1 1 dt' 

A, (x(t'), t') + B x (x(t') , t)ln-^pj I — 



and 

Av<?\x,t) = Av^\t' 2 ) exp 



(75) 



to 



A 2 (x(t'),t') + B 2 (x(t'),t')ln 



1 1 dt 7 
xW(t') J "i 7 " 



(76) 

These two equations (Eq. (J7Hj) and Eq. (|7fjjl ) give the two unknown functions 
At»W(a;, t) and Av^(x, t) which are related to the polarized singlet AS(x, t) 
and polarized gluon Ag(x,t) by Eq. (f47IJ) . Solving Eq. (|4f)]) for Av(x,t) we get 



and 



Av^(x,t) 



Av (2 \x,t) 



: b' n +b' 22 «) A ( t) _jk AE( t) 

2 s s 



" yil t^ 2 + g) A^,t) + ^AS(x,t) 
2 s s 



(77) 
(78) 



where s is given by Eq.(|5Uj). From Eq. (|77jl and Eq. (|75j) we get 



AE(x,t) 
and 



-b'u + b' 22 + s) 



2b' 



Av {1) (x,t) 



b'u + b' 22 - s) 



21 



26', 



21 



Ag(x,t) = Av (1 \x,t) + Av (2) (x,t). 
In the special case when x —>■ 0, 



A*/ 2) (x,t) (79) 
(80) 



AS(x,t) = -(-8 - v / 55)Aw (1) (x,t) + i(-8 + v / 55)Aw (2) (x, t) (81) 
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and 

Ag(x, t) = Av {1) (x, t) + Av {2) (x, t) . 
That is, at very low x approximation 



(82) 



AE(x,t) 



i(-8 
x exp 



to 



55)Av^(r[) 
A 1 {x{t'),t') + J Bi(a;(t , ),t)/n 
hh)Av^{r' 2 ) 



xW(t') 



x exp 

and 

Ag(x,t) = Av {1 \t[) exp 
+Av {2 \r 2 )exp 



to 



A 2 (x(t'),t') + B 2 (x(t'),t')ln 



xM{t') 



dt' 



dt' 



(83) 



to 



to 



|Ai(x(t'),*') + B x {x{t'),t)ln 
A 2 (x(t'),t') + B 2 (x(t'),t')ln 



1 



( 2 )(f) 



> dt' 

dt'" 



.(84) 



In Eq.flBSJ) and Eq.flBU), A?;W(t() and Av^ 2) (t^) are the initial conditions 
on the initial curve t = to- Since At>W and Av ^ are some combinations 
of AE and Ag as given by Eq.((77j) and Eq. (j78|) . therefore, from the input 
distribution for AE and Ag at Q 2 = Qq (t = t ) , we can get AuW(r{) and 
Ai/ 2 )^) simply by substituting x —> t[ and a; — * ^[20]. But in carrying 
out the integration we face one problem. Under the integral sign in Eq. (J83j) 

and Eq. (j84J) . we have the ratio —(see Eq. (l69H72*|) ) which is yet unknown. 

This ignorance forbids analytical forms for the quantities defined in these 
equations. However, at the initial scale t = t this ratio is known, because 
we know AE and Ag. We assume that the x and t dependent parts of 
Av^ and Av ( 2 > are factorisable and the x dependent part does not deviate 
significantly from their values at t = to i.e. 



Av^(x,t) Av^(r^ 



AvW(x,t) AvW(t[) 
where Af(t) is some unknown test function 



£A/(t) 



(85) 
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2.2 Analytical solution when t is very near the bound- 
ary t = t 

It is possible to obtain some analytical forms for the function AS and Ag 
in a region very close to the boundary t = t . Expanding Af(t) in a Taylor 
series about t = to and retaining only the first term we get 



Af(t) = Af(to) + (t-t )Af'(t) + 
» A/(* ) = A/ . 



(86) 



With this approximation the integration of Eq.(j75j) and Eq. (J76j) give the 
following analytical forms 



Av (1 \x,t) = Av^{r[) 



cxp 



and 



Av {2 \x,t) = Av {i \r' 2 ) exp 

where 
A 10 



B w In — In- 
to x 



B20 In — In- 
to x 



(87) 



1 



Bio 
B20 



2970 
1 

2970 
1 

2970 
1 

2970 



t ,(A 20 +B 20 ) 

- r'2^ x 1 Q < 2 ' } 

to J 



M T 2J 



(-2255 + 284V55) + 3(275 - 52V55)^Hr A /o 

M r i) 

(4400 - 404^) + 4(-880 + 119^55) ^^ A/ 



(2255 + 284V55) + 3(275 + 52V55) 



Av (1 \r[ 



Av^l 



T~2 



7V A /o 



A (1)/ /\ 

(4400 + 404^) - 4(880 + 119^) A ^ ){ A/ 

Av^>{T2) 

t[ = x + a {l) ln(t /t) 



(89) 
(90) 



and 



T ' 2 = x + a (2) ln(t /t) . (91) 

Using Eq.flBZD and Eq. (jBB)l in Eq.flESJ) and Eq. ipgjl . we get the analytical 
solutions for the polarized singlet and the gluon distribution which are valid 
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when x is very small and t is very close to the boundary t = to of perturbative 
evolution. Explicitly these are 



AE(:M) = ±(- 



8 - v / 55)Aw (1) (r()exp 



-Bio In — In- 
to x 



— I n-'\ ff ( „(1) J 



8 + v / 55)Aw (2) (r2)exp 



r, r t , 1 

-D20 ™ — ^ — 
t £ 



X 



and 

Ag(x,t) = Av (1) (t[) exp 
+Ai; (2) (t2) exp 



(92) 



-Bio In — In- 
to x 

n , t , 1 

-B20 *7l — 

t X 



-) r'p Tr x~ ( ^r ) 
to/ 

- r' 2 « (2) x ( Q < 2 » j . 

to J 



(93) 



3 Results and discussion 

In recent years there have been considerable growth of data fTJ 1221 
123 123 123 I21j on inclusive polarized deep inelastic scattering of leptons off 
nucleons. These data on the polarized asymmetry provide indirect way of 
extracting the polarized parton distribution functions that are essential to 
understand the spin decomposition of the nucleon. The distributions are 
determined from a global QCD analysis of polarized experimental data using 
the DGLAP evolution equations. Several such distributions are available in 
the literature El El El El HHj- We compare our results derived in the previous 
section with some of these exact solutions |5J El Ej of DGLAP equations that 
are available in the hep database as Fortran codes [2*7]. 

We first compare the polarized gluon distribution AG (= xAg) given 
by En. ipjl with the exact distribution AACOOjS] at Q 2 = 5GeV 2 . The 
analytical expression Eq. (j9l?|) . expected to be valid at very low x and low t has 
a free parameter A/ introduced through the Eq.([86|). Since A/ essentially 
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represents the ratio of some combination of singlet to the gluon polarized 
distribution at the initial scale, we assume that its value cannot be zero or 
negative. Starting from a low positive value, we vary this parameter till we 
get a reasonable agreement with the exact distributions. In Fig. 1(a) we show 
the graphs only for few representative values of A/o along with the exact 
AACOOLO0 distribution. To evolve our polarized gluon we also use the 
AACLO input at the initial scale Q 2 Q = lGeV 2 as given in the ref.jH]- From 
the figure(Fig.l(a)) we see that, while at x < 0.01, there is no appreciable 
difference between our distributions for different values of the parameter 
A/o, above x > 0.01, the difference is significant. We can see that for A/o ~ 
2.4, our prediction conforms reasonably well with the exact distribution. 
However, if we change the input distributions, then the value of the parameter 
A/o for which we get a reasonable agreement with the exact distribution 
also changes. In Fig.l(b-c) we compare our result with two other exact 
distributions: LSS01|S| and GRSVO10 at the same scale Q 2 = 5GeV 2 . For 
evolution of Ag we use the inputs also from the above two references. From 
the figure (Fig.l(b)(c)) we see that agreement between our prediction and the 
exact solution is good when A/ ~ 2.2 in the case of LSS input distribution 
and A/o ~ 2.0 in the case of GRSV input distribution. Thus we find that 
the value of A/o is sensitive to the input distributions and different input 
distributions give different A/ for a reasonable agreement with the exact 
ones. We use these values of A/ for the later part of our discussion to make 
the analysis self consistent. 

Now let us come to the gluon asymmetry Ap. Presently, the experimen- 
tal data available on the gluon asymmetry are limited and there are only 
three data sets from SMC^H, HERMES j23J EH 123 and COMPASS [213 I2E] 
obtained from the asymmetry measurement of high-p^ hadron production. 
These are shown in table 1. The bracket <> in the table means that the 
value is obtained at an average fraction of the nucleon momentum carried by 
the struck quark. The measurement is carried out at a momentum scale of 
the order of Q 2 = 2 ~ 5GeV 2 . We calculate the quantity AJ^, with AG given 
by Eq. (J93|) and the unpolarized gluon distribution taken from the solution 
of the LO DGLAP gluon evolution equation as obtained in |2H|. Unpolar- 
ized gluon distribution is evolved with MRSTLO[30J input distribution at 
Q 2 = lGeV 2 . We show our calculated values for A2 at Q 2 = 5GeV 2 along 
with the available data in Fig.2. The vertical error bar in the data represents 
the statistical error and the horizontal error bar, the standard deviation of 
the x g distribution at which the measurement is carried out. From the figure 
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Figure 1: Polarized gluon given by the analytic expression Eq. (j93|) compared with the exact (a)AAC00LO[5 , 
(b)LSSLO 6j and (c) GRSVLO[7] distributions at Q 2 = 5GeV 2 for different values of the free parameter 
A/ . 



Table 1: Recent Gluon asymmetry ^p data from COMPASS, SMC and 
HERMES 

COMPASS ^ = 0.06 ± 0.31(0 ± 0.06(syst) < x >= 0.13 
SMC ^ = -0.20±0.28(st) ±0.10(syst) < x >= 0.07 

HERMES 4r = °- 41 ± 0-018(st) ± O-OS(syst) 0.06 < x < 0.28 



we see that our predicted ^p reproduces the general trend of the values of ^p 
indicated by the data. We also see that the gluon asymmetry ^p invariably 
remains positive and conforms to the positivity constraint ^p < 1. However, 
this is true only at low momentum scale Q 2 < 5GeV 2 . As the momentum 
scale is increased, the condition -^p < 1 breaks down at certain values of x 
which we discuss in the next paragraph. 

At this point, few remarks about the positivity constraints are in or- 
der. The positivity condition has the origin in the probabilistic interpreta- 
tion of the parton densities. According to it, the magnitude of a polarized 
cross-section should be smaller than the corresponding unpolarized one, i.e. 
| Act | < a. In LO, probabilistic interpretation can be applied to parton dis- 
tributions and hence this condition (i.e. |A<r| < a) implies that the parton 
densities should satisfy the relation 

A/ (x,Q 2 )\ <f(x,Q 2 ). 

In determining the polarized distributions by different groups, the positivity 
condition is imposed along with other rules such as quark counting rulej2U 
E21 E2] to find the parameters of the initial distributions. In our study of the 
polarized DGLAP equations, since we are not making any fit except varying 
the parameter A/o, we see for what kinematic regions, our solutions conform 
to the positivity condition for a particular input distribution. Taking the 
LSSLO[H] input distribution, we calculate ^p for a large number of fixed Q 2 
and see for what value of x, if there is any, the condition ^p < 1 breaks 
down. The value of the free parameter A/o is chosen to be 2.2 (another 
graph for A/o = 2.0 also shown) as obtained from the previous analysis. In 
Fig. 3 we show the x — Q 2 range within which the condition Ap < \ holds 
good. We see that as Q 2 is increased, the condition Ap < 1 breaks down at 
progressively lower values of x. On the other hand, for low momentum scale 
Q 2 < 5GeV 2 it is always less than one. 
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Figure 2: Gluon asymmetry with two different inputs compared with data 
from SMC[221,COMPASSpni and HERMES .Values of A/ are chosen 
from figure [T]as discussed in the text. 
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Figure 3: x — Q 2 range where the positivity condition -£r < 1 holds good. In 
the region below the curve < l.The graphs are drawn with LSSjHj input 
in Eq. (|93jl for two different values of A/q. Unpolarized gluon G(x, t) is taken 
from [22]. 
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Next we consider the flavour singlet polarized quark distribution AS (x, t) 
given by Eq.()92)l and compare it with exact AACLO0 and LSSLO[U] distri- 
butions. Taking the inputs also from AACLOJH] and I.SSI.O (i respectively, 
we evolve AS(x, Q 2 ) using Eq. lj^j) at two different scales Q 2 = 2GeV 2 and 
Q 2 = 5GeV 2 and plot these as function of x along with the exact solutions 
in Fig. 4. The value of A/ is taken to be 2.2 for AAC input and 2.4 for LSS 
input respectively as determined earlier. From the graphs we see that the 
qualitative features of the LO exact distribution is well reproduced by the 
analytical form of AS (a;, Q 2 ) given by Eq. (J92J) . 

The value of AE(x, Q 2 ) and Ag(x,Q 2 ) also give us information about 
the contribution of the quark and the gluon to the total spin of the proton. 
The contribution of the singlet quark AS (x, Q 2 ) and the gluon distributions 
Ag(x, Q 2 ) to the total spin of the proton are given by their first moments: 



To calculate the first moment we have to integrate over the entire range 
of x from to 1. But in our derivation, we have an upper limit of x up 
to which we can calculate our functions AS(x, Q 2 ) and Ag(x,Q 2 ). This 
limit is set by the functions r[ and t' 2 above which they exceed the physical 
limit of unity. This upper limit is dependent on Q 2 . For Q 2 = 5GeV 2 , this 
upper limit is x max ~ 0.6. We calculate the first moments for our functions 
AS(s,<5 2 ) and Ag(x,Q 2 ) at Q 2 = 5GeV 2 and show them in table 2 along 
with the results obtained by GRSV[7J and AAC[5] LO parametrization. In 
table we have quoted only the LO result. The difference in their values for 
the spin content of the quark mainly originate from the small x behaviour of 
the anti-quark distributions and the fact that accurate experimental data are 
still not available at small x. The usual range of AS obtained by different 
parametrization groups lies between 0.1 ~ 0.3 and we see that our calculated 
AS is within this range. For the gluon polarization we have obtained a low 
value compared to the two parametrizations, presumably due to the cut-off 
imposed on x max in evaluating the integrals. 

For completeness, we also note that a parallel analysis is possible if one 
assumes an ad-hoc analytical form of the unknown function Af(t) in Eq. (j85J) 
instead of the parameter A/ (Eq.(jHHJ)) adjusted to the exact result. We 
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Figure 4: Polarized singlet distribution given by Eq. (j92j) as a function of x compared with the exact 
AACOOLO0 and LSSLO|E] distributions at two fixed Q 2 = 2GeV 2 and Q 2 = 5GeV 2 . Values of A/ are 
chosen as discussed in the text. 



Table 2: Quark and Gluon helicity distribution 



Q 2 (GeV 2 ) 


AAC 


GRSV 


ours 


AS 5 


0.18 


0.259(standard) 


0.201 






0.248(valence) 




Ag 5 


1.314 


0.684 (standard) 


0.31 






0.963 





refrain from doing such an analysis since Q 2 explored in the PDIS is limited 
near the boundary t rs t and effectively A/(i) ~ A/(t ) = A/ . 
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